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FRACTURE ANALYSIS OF SURFACE AND THROUGH CRACKS 


IN CYLINDRICAL PRESSURE VESSELS 

J. C. Newman, Jr. 

Langley Research Center 


SUMMARY 

A previously developed fracture criterion was applied to fracture data 
for surface- and through-cracked cylindrical pressure vessels to see how 
well the criterion can correlate fracture data. Fracture data from the lit- 
erature on surface cracks in aluminum alloy, steel, and epoxy vessels, and 
on through cracks in aluminum alloy, titanium alloy, steel, and brass ves- 
sels were analyzed by using the fracture criterion. The criterion corre- 
lated the failure stresses to within ±10 percent for either surface or 
through cracks over a wide range of crack size and vessel diameter. The 
fracture criterion was also found to correlate failure stresses to within 
±10 percent for flat plates (center-crack or double-edge-crack tension 
specimens) and cylindrical pressure vessels containing through cracks. 


INTRODUCTION 

Failures of many pressure vessels have been traced to surface cracks 
or to through cracks. These cracks initiate at structural discontinuities 
such as holes, material defects, or other abrupt changes in configuration 
and may propagate to failure under operating stress levels. To prevent such 
failures, the designer must be able to predict the effects of crack size on 
structural strength. Linear elastic fracture mechanics (LEFM), which uti- 
lizes the concept of the elastic stress-intensity factor, has been used to 
correlate fracture data and predict failure for cracked plates and struc- 
tural components when the crack-tip plastic deformations are constrained to 
small regions (plane-strain fracture, ref. 1). However, when plastic defor- 
mations near the crack tip are large (plastic zone greater than plate thick- 
ness) the elastic stress-intensity factor at failure varies with crack 

size and structural dimensions. (See refs. 2 to 4.) To account for the 
variation in with structural dimensions, the elastic-plastic stress- 

strain behavior at the crack tip must be considered. 

An equation which accounts for the effects of plastic deformation on 
fracture was derived and is presented in references 4 and 5. This equation 
relates to the elastic nominal failure stress and two material frac- 

ture parameters and is designated the two-parameter fracture criterion 
(TPFC). The TPFC has correlated fracture data for surface and through 
cracks, for different specimen types, and for a wide range of materials 
(refs. 4 to 7). 



The purpose of this paper is to apply the TPFC to surface- and through- 
cracked cylindrical shells subjected to internal pressure (fig. 1) to see 
how well the TPFC can correlate such fracture data. Fracture data for pres- 
surized cylinders made of various materials (aluminum alloy, titanium alloy, 
steel, brass, and epoxy) were taken from the literature. (See table I and 
refs. 8 to 17.) Some of the literature sources reported fracture data on 
flat plates (fig. 2) in addition to the data for cylinders. To determine 
whether or not the failure stresses in the cylinders may be predicted using 
fracture data for the flat plates, experimental and predicted failure 
stresses in the flat plates were compared . The predicted failure stresses 
were computed using the two material fracture parameters determined from the 
cylinder data. In order to apply the TPFC to fracture data, the elastic 
stress-intensity factors for these crack configurations must be known. 

The elastic stress-intensity factors for through cracks in pressurized 
cylinders have been obtained theoretically by Folias (ref. 18), and Erdogan 
and Kibler (ref. 19) for crack lengths less than about five times the square 
root of the product of vessel radius and thickness. In this paper an empir- 
ical equation giving elastic stress-intensity factors for through cracks has 
been obtained using the numerical results from reference 19 and some experi- 
mental fracture tests on brass vessels from reference 16. The resulting 
equation applies over a range of crack lengths about twice as large as those 
considered in references 18 and 19. 

The elastic stress-intensity factors for surface cracks in pressurized 
cylinders have not been obtained theoretically, but some experimental stress- 
intensity factors have been determined from surface cracks in brittle epoxy 
vessels (ref. 17). In this paper an empirical equation giving elastic 
stress-intensity factors for surface cracks has been developed using the 
results from references 4 and 19. The results from this equation are com- 
pared with the experimental stress-intensity factors from reference 17. 


SYMBOLS 

a initial depth of surface crack, m 

c Initial length of surface or through crack (see figs. 1 and 2), m 

F complete boundary-correction factor on the stress intensity 

f^ shell-curvature correction factor for a through crack 

f shell-curvature correction factor for a surface crack 

s 

H height of uniformly stressed specimens, m 

Kp fracture toughness , N/m^/^ 

Kj elastic stress-intensity factor, 

Kjg elastic stress-intensity factor at failure, 
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Mg combined front-face and back-face correction on the stress- 

intensity factor for the surface crack 

M-| front-face correction on the stress-intensity factor for the 

surface crack 


m fracture-toughness parameter 


p Internal pressure, Pa 

Q elastic surface-crack shape factor 

R internal radius of cylindrical pressure vessel, m 

Sg gross section stress at failure. Pa 

nominal (net section) stress at failure. Pa 

Sy nominal stress required to produce a plastic hinge on net section 

f’of’ center-crack and double-edge-crack tension speci- 
mens and = 1.15a^ for pressurized cylinders). Pa 

T temperature, K 

t shell thickness, m 

W specimen width, m 


through-crack shell parameter, 
surface-crack shell parameter. 



Poisson's ratio 

ultimate tensile strength (uniaxial). Pa 
yield stress (uniaxial). Pa 
ratio of to 


TWO-PARAMETER FRACTURE CRITERION 

The two-parameter fracture criterion (TPFC) was developed and success- 
fully applied to plane fracture specimens containing either surface or 
through cracks in metallic materials (refs. 4 to 7). The TPFC accounts for 
the effects of plastic deformation on fracture properties. The equation is 
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( 1 ) 


K 


F 





where is the elastic stress-intensity factor at failure, is the 

nominal tnet-section) failure stress, is the nominal stress required 

to produce a plastic hinge on the net section using the ultimate tensile 
strength, and K„ and m are the two material fracture parameters. The 
fracture parameters K„ and m are assumed to be constant for a given com- 
bination of material, thickness, temperature, and rate of loading. Three 
conditions are necessary to obtain fracture parameters that are representa- 
tive for a given material and test temperature: (1) the nominal failure 

stress must be less than a , (2) the fracture data must all be from the 
same specimen thickness, and (3) the test data must encompass a wide range 
of crack lengths or specimen sizes. Reference M shows how the fracture 
parameters are determined by a least— squares procedure for a given set of 
fracture data. To define the complete fracture behavior for a material, Kp 
and m must be determined as functions of thickness, temperature, and load 
rate . 


If m equals zero in equation (1), Kp equals the elastic stress- 
intensity factor at failure and the equation applies to low-toughness (low 
Kp) materials (plane-strain fracture). However, if m equals unity, the 
equation applies to extremely ductile or high-toughness (high Kp) materi- 
als. Thus, the fracture parameters, Kp and m, jointly describe the crack 
sensitivity of the material. 

The denominator in equation (1) reflects the influence of the nominal 
f'gilure stress on fracture toughness. The variation of the denominator with 
nominal stress for a typical material is shown in figure 3. When the nomi- 
nal stress is less than the uniaxial yield stress , the function <t) 

(ratio of K, to Kp) is a linear function of nominal stress. However, 
when the nominal failure stress is greater than the yield stress, the func- 
tion (j) becomes nonlinear and is dependent upon the stress-strain curve of 
the material and the- state of stress in the crack-tip region (refs. 4 
and 6). For thin materials, where the state of stress in the crack-tip 
region is biaxial, the expected behavior is estimated by the dash-dot curve. 
An equation approximating the dash-dot curve is given in reference 6 by 


<t> 




(cj,. < S, 




(2) 


and is shown in figure 3 as the dashed curve. For thick materials, where 
the state of stress in the crack-tip region is triaxial , the fracture behav- 
ior for S ^ is expected to lie closer to the solid line. The solid 
vertical line truncates the nominal stress at In this paper, equa- 
tion (2) was used whenever was greater than 
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In order to apply equation ( 1 ) to surface and through cracks In pres- 
surized cylinders, the nominal stress required to fail the uncracked vessel 
and the elastic stress-intensity factor Kjg for these configurations 
must be determined. For pressurized cylinders, was 1.15 times the ulti- 

mate tensile strength u . It was determined by calculating the nominal (hoop) 
stress required to satisfy the Mises yield criterion (where cr^ was replaced 
by cr^) assuming a 2:1 biaxial stress ratio. The elastic stress-intensity 
factor equations for these configurations are presented in the next section. 


ELASTIC STRESS-INTENSITY FACTORS 

The form of the elastic stress distribution near a crack tip that con- 
tains the stress-intensity factor and the square-root singularity is 

well known (ref. 2). (The determination of is the basis for linear 

elastic fracture mechanics.) The stress-intensity factor is a function of 
load, structural configuration, and the size, shape, and location of the 
crack. In general, the elastic stress-intensity factor can be expressed as 

Kfe = Sn^ F (3) 

for any Mode I crack configuration where S is the nominal stress, c is 
the initial crack length (defined in figs. T and 2), and F is the boundary- 
correction factor. The boundary-correction factor accounts for the Influence 
of various boundaries and crack shape on stress intensity. The following 
sections give the nominal stress equation and the boundary-correction factor 
equations for the surface- and through-cracked cylindrical pressure vessel . 


Through Crack in a Thin Pressurized Cylinder 


For the through crack (axial) in a cylindrical shell subjected to inter- 
nal pressure (fig. 1(a)), the elastic stress-intensity factor at failure is 
given by equation (3) where the nominal stress is 


and 



(4) 


F = ft = (l + 0.52^^ + 1.29'^^ - 0.074^^3 (5) 

for 0 ^ 10 where = c/^Rt. Equation (5) accounts for the effects 

of shell curvature (refs. 18 and 19) on stress intensity. Poisson's ratio 
was assumed to be 1/3- The details on the development of equation (5) are 
given in the appendix. 


Surface Crack in a Thin Pressurized Cylinder 

For the internal or external surface crack (axial) in a cylindrical 
shell subjected to internal pressure (fig. 1(b)), the elastic stress- 
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intensity factor at failure is also given by equation (3), where the nominal 
stress is given by equation (4) and 



The square-root term converts the through-crack expression to that for a 
surface crack, is the combined front-face and back-face correction fac- 

tor, and fg is the shell-curvature correction factor for a surface crack. 
The elastic shape factor Q was given in reference 20 as the square of the 
elliptic integral of the second kind. An expression was chosen in refer- 
ence 4 as a simple approximation for Q and is given by 


Q 



Q = 




( 7 ) 


The expression for (ref. 4) is given by 



( 8 ) 


where q was determined empirically as 

, = 2 * 


(9) 


The term is the front-face correction, and the a/t term is the back- 

face correction. The expression for is given by 


M. = 1.13 - 0.l(- 
I \c/ 




0.02 - 4 1-0 
c 


( 10 ) 




c\ 

/a \ 

= \l 

- 1 + 0.03 - 



lav 

a/ 

Vc / 

J 


6 


is given by 


The shell-curvature correction factor for a surface crack f 

s 

fg = (^1 + 0 . 52 X 3 + 1.29Xg2 - O.OTllXgSj (11) 

C 3. 

for 0 ^ ^g ^ 10 where Xg = - — Again, Poisson's ratio was assumed 

/Rt ^ 

to be ^/ 3 ^ The form of x^ was obtained by assuming that the surface 
crack could be replaced by an "equivalent" through crack of equal area. 

As a/t approaches unity, Xg approaches X^ and equation ( 6 ) reduces 
to equation (5). In the appendix, equation ( 6 ) is compared with some 
experimentally-determined correction factors for a brittle epoxy. 


Through Crack in Axially Loaded Flat Plates 

C enter-crack tension .- For the center-crack tension specimen 
(fig. 2 (a)), the elastic stress-intensity factor is given by equation ( 3 ) 
where 


S 


n 



( 12 ) 


and 


(13) 


The secant term is the finite-width cor- 
rection on stress-intensity factor and was obtained from reference 1 . 

Dpuble-gd ge-crack tension .- For the double-edge-crack tension specimen 
(fig. 2 (b)), the elastic stress-intensity factor is given by equation ( 3 ), 
where is given by equation (12), and F was obtained by a boundary- 

collocation analysis of a configuration with H/W = 0.625. (Fracture data 
analyzed from ref. 15 used this particular configuration.) The stress S 
was assumed to be uniformly applied. An equation fit to the collocation ® 
results gave 


F = 


. 2C\ / TTC 

1 - — I l/sec — 

W / H W 


for 0 £ 2c/ W < 1.0 and H/W 2 2, 


F 


1.13 1 



(14) 


for 0 ^ 2c/W ^ 0.4 and H/W = 0.625. (Stress-intensity factors for 
H/W =1,3, and ” are given in ref. 2.) 
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ANALYSIS OF FRACTURE DATA 


Fracture data on pressurized cylinders made of various materials and 
containing either surface or through cracks were taken from the literature 
and were analyzed using the TPFC (eqs. (1) and (2)). The fracture param- 
eters, Kp and m, for a given material, thickness, and test temperature 
were determined from the fracture data using a best-fit procedure described 
in reference 4. These values of Kp and m were then used to calculate 
failure stresses for the same fracture tests to see how well the TPFC cor- 
related the failure stresses. Where fracture data on flat plates (center- 
crack tension or edge-crack tension) were available, experimental and pre- 
dicted failure stresses for the flat plates were compared. The failure 
stresses were computed using the two material fracture parameters determined 
from the cylinder data. The failure stresses were calculated by substitut- 
ing equation (3) into equations (1) and (2), and were given by 


S 


n 


K 


F 


mKp 



(Sn 1 <^y) (15) 


and 


+ 2YSy - mY 






( 16 ) 


where 


2S^^ F 

Table I summarizes the materials, thicknesses, test temperatures, and crack 
configurations that were analyzed. Figures 4 to 15 show the correlation 
obtained for each set of fracture data. Table I also indicates which frac- 
ture data were obtained using sharp saw-cut slits instead of fatigue cracks. 


Through Cracks 

Pressurized cylinders .- Figures 4 to 7 show fracture data on pressur-^ 
ized cylinders containing through cracks. The figures show the nominal fail- 
ure stress normalized to S (1.15c ) plotted against half-length of crack 
c. The symbols show the fracture data, and the curves indicate the best fit 
of the TPFC using the values of Kp and m determined from these data. A 
knee occurs in all curves when the nominal stress is equal to the yield 
stress of the material (transition from eq. (15) to eq . (16)), but some are 
hardly perceptible. The calculated failure stresses were generally within 
±10 percent of the experimental failure stresses. 
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Pressurized cyl inders and flat plates .- Figures 8 to 13 show fracture 
data on pressurized cylinders and flat plates. In each figure the results 
are for cylinders and flat plates made of the same material and thickness, 
and tested at the same temperature. Both the cylinders and flat plates con- 
tained through cracks. The figures show nominal failure stress normalized 
to Sy plotted against half-length of crack c. For the flat plates 
^u ” % cylinders = 1.15^ The fracture parameters, Kp, 

and m, were determined from an analysis of the fracture data on the pres- 
surized cylinders. The fracture data on the flat plates were insufficient 
to obtain the two parameters because only one specimen size and crack length 
was tested . The solid curve or curves show the calculations from the TPFC 
for various vessel radii. The agreement was considered good. The dashed 
curve on each figure shows the predicted results for flat plates (R = “) 
using the values of Kp and m determined from the cylinders. The circu- 
lar symbols show the experimental results for the flat plates. The pre- 
dicted results were within ±10 percent of the experimental failure stresses. 


Surface Cracks 

Pierce (ref. 11) conducted cryogenic fracture tests on surface cracks 
in 2014-T6 material. Figure shows normalized to S (1.15cJ ) 

plotted against cF^, where c is the half-length of crack'^and F is 
the boundary-correction factor (eq. (6)). The TPFC indicates that this 
type of plot against cF^) gives a single curve for various a/c 

and a/t ratios. The open and solid circular symbols denote experimental 
data on externally or internally located surface cracks, respectively. 

The surface-crack data included variations in crack shape [o.07 £ - 1 0.9) 

f a \ \ ° / 

and crack size 10.36 S - ^ 0.98). The curve shows the calculations from the 

TPFC using the values of Kp and m determined from these data. The cal- 
culated failure stresses were within ±10 percent of the experimental failure 
stresses . 


Kiefner , Maxey , Fiber, and Duffy (ref. 1^4 ) conducted surface-crack frac- 
ture tests at room temperature on steel pressure vessels. Figure 15 shows 
the nominal failure stresses normalized to plotted against cF^ . The 

surface-crack data included variations in crack shape ( 0.02 < - = 0.1 4 and 

/ a \ \ ° / 

crack size (0.38 1 - ^ O. 82 I. The symbols denote the fracture data and the 

the curve shows the calculations from the TPFC. The agreement was consid- 
ered good . 


CONCLUDING REMARKS 

A two-parameter fracture criterion that relates the elastic stress- 
intensity at failure, the elastic nominal failure stress, and two material 
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parameters was used to analyze fracture data on surface- and through-cracked 
cylindrical pressure vessels. Fracture data from the literature on steel, 
titanium alloy, aluminum alloy, brass, and epoxy vessels tested at either 
room or cryogenic temperature were analyzed. The two-parameter fracture cri- 
terion correlated the data well (generally within ±10 percent of the experi- 
mental failure stresses) for a broad range of materials, including some that 
were extremely ductile. The fracture criterion was also found to correlate 
fracture data from flat plates and cylindrical pressure vessels within 
±10 percent for the same material, thickness, and test temperature. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
September 27 , 1976 
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APPENDIX 


DEVELOPMENT AND VERIFICATION OF BOUNDARY CORRECTIONS ON STRESS-INTENSITY 
FACTORS FOR THIN PRESSURIZED CYLINDERS 
Through Cracks 

Follas (ref. 18) and Erdogan and Kibler (ref. 19) have obtained the 
elastic stress-intensity factors for a longitudinal (axial) through crack in 
a pressurized cylinder (fig. 1(a)) for ^ 4.5. Figure 16 shows the shell- 
curvature correction (or boundary-correction factor) on stress intensity as 
a function of x^ . The symbols are numerical values obtained from refer- 
ence 19 for V = 1/3. The shell-curvature correction factor F was deter- 
mined by 



where F-j^, is the contribution due to the membrane solution (or normal 
forces) and Fg is the contribution due to bending. Reference 21 has shown 
that including one-half of the bending contribution was necessary to corre- 
late crack-growth rates from pressurized cylinders and flat plates. The 
bending term contributed only about 10 percent to the total correction 
factor . 

The solid curve in figure 16 is an equation chosen herein (eq. (5)) to 
fit the numerical values. The solid curve was within ±5 percent of the 
numerical values. Reference 21 has shown that the effect of Poisson's 
ratio ( V = 0 to 1/2) on the curvature correction was less than 5 percent 
from the values given in figure 16. Therefore, equation (5) was assumed to 
apply for any material. 

For X^ >4.5, an experimental technique was used herein to verify the 
applicability of equation (5) for calculating the shell-curvature correction 
factors. Figure 17 shows the shell-curvature correction factor plotted 
against X^. . The symbols show the experimentally derived correction factors 
from through-crack fracture data on brass cylinders (ref. 16) using the two- 
parameter fracture criterion. The correction factors are given by 



The fracture parameters Kp and m were determined from an analysis of 
the fracture data with X^ ^ 4 . The crack length c and the corresponding 
failure stresses were obtained from reference 16. The solid curve 

shows the correction factors calculated from equation (5), and these factors 
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are in good agreement with the experimental data. Therefore, equation (5) 
was assumed to apply for any material with £ 10. 


Surface Cracks 


Derby (ref. 17) has experimentally determined stress-intensity correc- 
tion factors for surface cracks in pressurized cylinders made of a brittle 
epoxy. The critical elastic stress-intensity factor for this material 
(Kj = 1.02 MN/m^^^) was obtained from four-point notch bend fracture tests 
(ref. 17). This material was brittle (m = 0). The experimental correction 
factors for the surface cracks were obtained from equation (3) as 


%e 1.02 



( 20 ) 


where S and c were obtained from the fracture tests on the epoxy ves- 
sels. Table II shows a comparison between the experimentally determined 
boundary correction factors and those calculated from equation (6). The 
theoretical correction factors were within ±10 percent of the experimental 
values . 
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TABLE I.- MATERIALS, THICKNESSES, TEST TEMPERATURES, AND CRACK CONFIGURATIONS 


I 


Material 

. Thickness , 
mm 

1 m 

Teroper- 
1 ature , 

I ' 

i K 

! Pressurized 

L 

cylinders 

Flat plates, 
through 
crack 

Reference 

Surface crack 

Through crack 

Aluminum alloy: 







7075-T6 

0.4 to 0.5 

: ^RT 




8 

2024-T3 

0.3 to 0.4 

PT 




8 

2014-T6 

1.5 

RT 




9 

2014-T6 

1.5 

20 


X 

X 

10 

2014-T6 

1 .5 

77 


X 

X 

10 

2014-T6 

1.5 

77 

X 



11 

Titanium alloy: 







Ti-5Al-2.5Sn (ELI) 

0.5 

20 


X 

X 

10 

Ti-5Al-2.5Sn (ELI) 

.5 

77 


X 

X 

12 

, Steel : 







AISI 301 

0 . 6 

RT 


X 


13 

X-52 

9.5 

RT 


^x 


14 

X-52 

9.5 

RT 




14 

Hot-rolled° 

6.4 

77 


X 

X 


Brass 

0.025 

RT 


^x 


16 

Brittle epoxy 

15 

RT 

X 



17 


Room temperature. 

*^Sharp saw-cut slits instead of fatigue cracks. 
°0.25C, 0.02Si, 0.85Mn. 




TABLE II.- THEORETICALLY AND EXPERIMENTALLY DETERMINED 
BOUNDARY-CORRECTION FACTORS EOR SURFACE CRACKS IN 
PRESSURIZED CYLINDERS MADE OF A BRITTLE EPOXY 

[R = 68.3 tnm and t = 15 mm (from ref. 17)] 


Theoretical F 

Theoretical F 
Experimental F 

0.932 

0.97 

.863 

.98 

.855 

1 .01 

.823 

.97 

.841 

1 .04 

.839 

1 .01 

.819 

1 .04 

.696 

1 .02 

.689 

.97 

.689 

1.03 

.677 

1 .09 
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Figure 3-- Typical 
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Figure 7.- Nominal failure stresses for through crac 
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Figure 8.- Nominal failure stresses for through cracks in flat plates and pressurized cylinders 

made of 2014-T6 aluminum alloy. 
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Figure 10.- Nominal failure stresses for through cracks in flat plates and pressurized cylind 

made of Ti-5Al-2.5Sn (ELI) titanium alloy. 
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Figure 17.- Experimental and theoretical shell-curvature correction factors 
for an axial through crack in a pressurized cylinder. 
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